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Abstract 

A tent-shaped structure with a flexible photovoltaic blanket acting as a catenary 
collector is presented. The shadow cast by one side of the collector on the other side 
producing a self shading effect is analysed. The direct beam, the diffuse and the albedo 
radiation on the collector are determined. An example is given for the insolation on the 
collector operating on Mars surface for the location of Viking Lander 1 (VL1). 
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1. INTRODUCTION 


Mission to Martian surface will require electric power. A power supply that requires 
little installation time* being light weight and stowed in a small volume can be accomplished 
by a photovoltaic (PV) array. A tent-shaped structure with a flexible PV blanket for solar 
power generation is proposed in [1]* Fig. 1. The array is designed with a self-deploying 

mechanism using a pressurized gas expansion. The structural design for the array uses a 

combination of cables* beams and columns to support and deploy the PV blanket. The array 
is stowed with the blanket either folded or rolled. The main contribution to the stress in the 
structure is due to the tension in the cable which supports the PV blanket. The shape of the 
PV blanket is determined by an optimization between reduction in the cable tension and 
increase in blanket area. Under the force of gravity a cable carrying a uniformly 

distributed load will take the shape of a catenary curve, f c (y), with respect to the Y-Z 

plane, Fig. 2. 

According to [2] and with some manipulations, f c (y) is given by: 


f c (y) 





2 ) 


l 


(D 


The catenary constant k can be determined using the condition f c (0) = H and solving 

iteratively. However, when the blanket is fairly taut, the load may be assumed uniformly 
distributed along the y-axis and the catenary curve may be approximated by a parabola [2], 
i.e.,: 
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fc(y) = §2 (D - y ) 2 


( 2 ) 


which simplifies the calculation. 

Because of the shape of a catenary-tent-collector, a self-shading effect occurs on one 
of its sides (side B in Fig. 2). In this article we analyze the shadow shape and area. Based 
on these results, the beam insolation on the collector is calculated. We also determine the 
diffuse and the albedo insolation on the collector. An example for the planet Mars is given 
for the location of Viking Lander 1. Results for the parabolic approximation is given in 
Appendix B. 


2. SHADOW CALCULATION 

The catenary-tent -collector is self shading. The size of the shadow and the side which 
is shaded depend on the sun position. In general, both sides of the collector will be 

alternately shaded in a given day if at sunrise 

“>sr < -1c - 90° (3) 

where 

~i c - the collector azimuth 
nf sr - the sun azimuth at sunrise. 

The azimuth angles are measured from true south positively in a clockwise direction. In days 
when eq. (3) is not satisfied, only one side of the collector is shaded all the times. In this 

paper we analyze the shape and size of the shadow cast on a catenary-tent-collector facing 

the south-north direction. The results can be generalized for an arbitrary oriented tent by 
replacing the sun azimuth angle n s by the difference between the solar and the collector 


3 


azimuth angles, i .e. , 


^ = ~is - ~*c ( 4 ) 

With reference to Fig. 3, the shadow cast by line ON can be divided into three distinct 
cases: 


(i) 

Pi < L 

. Py 

< 

D 


(ii) 

Pi < L , 

Py 

> 

D 

(5) 

(iii) 

Px > l 

. Py 

< 

D - 



It will be shown later that for the case P x > L, P y > D the shadow takes the shape either 
of case (ii) or case ( i i i ) . 

The components P z and P y of the line OP was derived in [3,4] and are: 

P z = H cos5 sincu/sina (6) 

and 

P y - H($in0 cos£ cost*; - cos0 sin£)/sina (7) 

where 

sina = sin0 sinS + cos0 cosS cosw (8) 

and 

8 - solar declination angle 
0 - local latitude 
a - sun elevation angle 
u - solar hour angle 
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Case «): Pr < L , Pw < D 


In this case the solar rays penetrate the collector surface and the shadow on the 
collector takes the shape NTEEjM as shown in Fig. 4. Point E is a point where the ray 
penetrates the collector, and point T is a point where the ray is tangent to the collector. In 
order to calculate the shaded area, the components of the points E and T must be 
determined. Since the catenary function f c (y) is hyperbolic, the solution for the component 
y E is obtained numerically from the-solution of the following equation (see Appendix A): 

H c ° sh ( e Tr“) + pf 1 1 + H (9) 

The component is then given by (see Fig. 4) 


Point T is in y-z plane, i.e., Xr = 0 and y r is obtained from the solution of (see Appendix 

A): 

yr = D - k sinh _1 (H/P y ) (11) 

It is easier to obtain the shaded area A S h from the "streched out’ collector, Fig. 4(b), i.e., 
Ash = L N E y - 0.5 Xr T E y (12) 


The unshaded area A ush is given by: 

A ush = A - Ash = L(N N n - N E*) + 0.5 X£ T E y 
where A is the area of one side of the collector, i.e., 


(13) 
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(14) 


A 


L N N 


n 


Defining a relative shaded area by 4, we obtain 
4 = (2 L - x E T - E y )/2L N~N„ 

The length of an arc between two points y A and yz is given by: 
yPy2 = f 2 -/T + [f c ’(y>] 2 dy 

J Vi 


(15) 


(16) 


where f c *(y) is the derivative of the function f c (y). The length of the arcs in eq. (15) thus 
become: 

N~N n = [ D Vl + Ifc’(y)] dy - k sinh(D/k) (17) 

Jo 


N Ey = J* E VT + [fc’(y)] 2 dy = k sinh(D/k) + sinh ^^ - jj (18) 


T Ej, = j E N /l + f c ’(y) 2 dy = kj^sinh ^ k + sinh^P yj )] (19) 


Cose (?7): P T < L , P v > D 

This case is shown in Fig. 5. The array touches the collector at the foot, point E, and 
at point D. The coordinates of point E are: 

Ye = D (20) 

As derived in Appendix A, 
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( 21 ) 


Xf = jr (D - yv) 


where y^ is determined from the solution of (see Appendix A): 


k [ c ° sh (D__JLu) - l] + » y„ 


DH 

P« 


( 22 ) 


The coordinates of point T are the same as for case (i), i.e., Xj- = 0 and yr is given by 
eq. (11). The shaded and unshaded areas, Fig. 5(b), are, respectively: 


A sh = A - 0.5 X£- T N„ 


A ush = 0.5 Xf T N„ 


(23) 

(24) 


The relative shaded area is obtained as: 


» _ , X£- T N n 

t - 1 - — 

2 L N N n 


(25) 


where the arc length T N n is given by: 


t"n„ = I ° y/1 + (fc’(y)) 2 dy = k sinh(V- T ) 


(26) 


and N N „ is given in eq. (17). 


Case (in): Pj > L * P v < D 

The solar ray may either intersect the collector or may fall outside the collector. The 
case of intersection was already analysed in case (i). To determine whether the solar ray 
intersects the collector* eqs. (9) and (10) are solved. 
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If xe < L than there is an intersection. If Xe > L than there is no intersection, see 


Fig. 6. The coordinates of point E are: 
Xf = L 


(27) 


and ye is determined from the solution of (see Appendix A): 


cosh 


D - y c + L 


V 


Px\ 


/ 


- coshl 


(V s ). 


L H 


(28) 


The coordinates of point T are = 0 and is given by eq. (11). The shaded and 
unshaded areas, Fig. 6, are respectfveiy: 


A S h = L N E y - 0.5 L T E y 


Aush = L(N N n * N E y + 0.5 T E„) 


(29) 

(30) 


The relative shaded area is: 

4 = (2 N~E y - TE y )/2 n"n„ (31) 

/■s / — s 

where N N„, N E y and T E y are given in eqs. (17)-(19). 

The case > L , P y > D corresponds either to case (ii) or to case (iii) depending 

on the ratio P*/P y of the ray OP (see Fig. 7). 

For Pj/P y < L/D the solution is as for case (ii) and 

for Pj/P y > L/D the solution is as for case (iii) 

For a south-north oriented collector, the shadow as a function of time is symmetrical 
with the solar noon. For the after noon, the shadow calculations are with respect to line 
O'M (see Fig. 3). 
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3. SOLAR RADIATION CALCULATION 


With the results developed in the previous section, the beam irradiance in W or W/m 2 
and the beam insolation in Whr/day or Whr/m 2 -day on the catenary-tent-collector can be 
determined. The diffuse and albedo components will be added to the beam to get the global 
irradiance and insolation. A north-south facing catenary-tent-collector will be considered. A 
generalization to any arbitrary oriented tent may be obtained using eq. (4). 

Beam Irradiance 

The beam irradiance on both sides of the tent depends on the self-shading condition of 
eq. (3). In calculating the irradiance on an unshaded side we resort to Fig. 8. The beam 
irradiance P, in watts (W), on an unshaded side is given by: 

dP* = G b cos 0 dA V (32) 

where is the beam irradiance in W/m 2 normal to the solar rays, and 6 is the angle 
between the solar rays and the normal to the surface, cos# and dA are given by 


xos0 = cosj3(y) sina + sinj3(y) cosa cosnr s (33) 

and 

dA = L dS = L y/l + [fc'(y)] 2 dy (34) 

where dS is the unit length of the collector (Fig. 8) and 

sin/3(y) = |fV(y)l/VT + [fc'(y)l 2 (35) 

cos/3 (y) = 1/v/l + [f c ’(y)] 2 (36) 
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Therefore, eq. (32) reduces to: 


dP b - G b L(sina + |f c '(y)| cosa cos^)dy (37) 

Integrating eq. (37) in the interval [0,D], we obtain the beam irradiance, in W, on the 
unshaded side of the collector, i.e., 

- G*> L(D sina + H cosa cos^) 38) 

Using the angle e as defined in Fig. 8 , we obtain: 

Pfc = G* L(H 2 + D 2 ) 1/2 (cose sina + sine cosa cosnr s ) 39) 

i.e., the beam irradiance on a catenary collector is equivalent to that of a flat plate 
collector MNM S N S . This conclusion is also valid for any non-flat shape collector. 


The beam irradiance on a collector that is partially shaded can not be obtained by 
multiplying the incoming beam irradiance with the factor ( 1 - 4 ) as for a flat plate collector 
since the beam irradiance on the collector is not uniform along the y axis of the collector, 
i.e., the angle 0 varies with y. Using eqs. (33), (35) and (36), the variation of the beam 
irradiance, in W/m 2 , along the collector is: 

Gi,(y) = 7 = ===== [sina + cosa cos-> s |f c ’(y)| ] 

VI + [fc (y)) 2 


and with eq. ( 1 ) we obtain: 


G/>(y) 



sina + cosa cos't s |sinh 



(40) 


With reference to Fig. 9(a), the beam irradiance on the partially shaded side of the 
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collector for case (i) is given by: 


Pf = P Gi(y) dAi + f G fc (y) dA 2 

Jy T •'Ve 


where G*(y) is given by eq. (40), and 


dA ! = x dS = x coshl —r~^ ) dy 


m 


( 41 ) 


dA 2 = L coshl - r -*- dy 


m 


(42) 


where x and y are related by (Fig. 9(a)) 


= (*^Vr) (y - yr) 


(43) 


Substituting these results, we obtain the beam irradiance on the partially shaded side of the 
collector: 


P£ h = Gi, { f [sina + cosa cos~r s | 


-*M fe) < y * »> < y 


f D r 

+ L / sina + cosa cosnr s ! 
■'2'e L 


i 


(45) 


Similarly, the beam irradiance for case (ii) (Fig. 9b) is: 


Pft h = G b j |sina + cosa cos*T s jsinh^^^ j|| (y - Yr) dy (46) 


and the beam irradiance for case (iii) (Fig. 9c) is: 
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= G b I J |sinu + cosa C05"i s |sinh| ^J|| ( yc - y T ) (y “ y r) 


/'[• 

r L 


+ L / j sina + cosa cosnr s | 
've 


sinh (nr)|] dy } 


( 47 ) 


In winter, side A (Fig. 2) remains always unshaded, therefore the beam irradiance P/m 
on this side is given by eq. (38) or eq. (39). Side B may be partially shaded and the beam 
irradiance P*£ is given by eqs. (45), (46) or (47) where the collector azimuth is 180°, 
therefore, ~i s is replaced by (nf s -180°). The total beam irradiance on the 
catenary-tent-collector in winter is: 


?b = p tA + n h B (**) 

In summer, early in the morning and late in the afternoon, side A is partially shaded, 
therefore, the beam irradiance P b A is given by eqs. (45), (46) or (47). Side B is unshaded, 
and the beam irradiance P b o is given by eq. (38) with the azimuth nr s -180°, i.e., 

P b B = G b L[D sina + H cosa cos(-? 5 -18G°)] (49) 

thus the beam irradiance on both sides of the tent collector for early morning and late 
afternoon hours is: 


P b - P t,A + P bB 


(SO) 


For other hours of the day, side B will be partially shaded and eq. (48) applies. 
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Diffuse Irradiance 


The diffuse irradiance on an incremental area dA is given by 

dP rf = G dh dA F rf( ,- S (51) 

where 

Gdh ~ diffuse irradiance on a horizontal surface 
F dA -s - view factor of dA with respect to the sky 

For sufficiently small dS, the view, factor for dA is given by: 

F^-s = 0.511 + cos£3(y)] (52) 


Using eqs. (34), (36) and (52) for eq. (51) we obtain: 


dP rf = 0.5 G rfh L(Vi + [f c ’(y)] 2 + Ddy 


Integrating for y in the interval [0,D] we get the diffuse irradiance on one side of the tent: 


P rf = 0.5 G dh L(N N„ + D) = 0.5 G dh L[k sinh(D/k) + D] 


(53) 


Since the diffuse irradiance is independent of the orientation of the collector (for isotropic 
skies), eq. (53) applies also to the other side of the tent. 

The view factor of a catenary collectors can now be calculated by defining: 


P</ = G dh A F^-s 

using eqs. (17) and (53) one obtains: 


F a - S = 0.5 


1 + 


D 


k sinh(D/k) 


(54) 


(55) 
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Albedo 


The albedo can be determined by using the expression: 

P„, = at G h A F A -c (56) 


where al - the albedo factor 

G* - global irradiance on a horizontal surface 
F A-c - view factor of area A with respect to ground 

The view factor of area A with respect to the sky, F^-s, and with respect to the ground, 
F A - a , are related by [5]: 

F A -s + F A -c = 1 (57) 

obtaining with eq. (55): 

F<4 ~ c = °- 5 [* ■ k sinh(D/k)] 

Using eq. (14), (17), (56) and (58) 
as: 

P a , = 0.5 al G h L[k sinh(D/k) 

Since the albedo irradiance is independent on the orientation of the collector eq. (59) applies 
also to other side of the tent. 


(58) 

we obtain the albedo irradiance on one side of the tent 
- D] (59) 


4. Example 

The example refers to a south-north facing catenary-tent -collector deployed on the 
Mars surface [6] at the locations of Viking Lander VLI (Latitude - 22.3°N % Longitude - 
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47.9°W) and in an autumn day L s = 200° [6], For this time of the year eq. (3) does not 
apply and, therefore, only side B will be partially shaded during the day. The dimensions of 
the catenary collector are: D = 3m, H = 2m and L = 1.5m. Since f c (0) = H, this results 
in k = 2.53 and the catenary equation is given by: 

f c (y) = 2 . 53 [ cosh C 1 - 1 86 - 0.395y) - 1] (60) 

Figure 10(a) shows the shadow shapes on side B of the catenary collector for 0 = 22.3°N 
for different hours of the day. The shadow calculations are based on section 2. For 
comparison, the shadow shapes for 0 - 22 ° are shown in Fig. 10(b). As expected, the 
self -shading effect is more pronounced for higher latitudes. It is interesting to note that for 
0=22. 3°N, the shadow effect is quite small during the noon hours. In summer, the shadow 
will be even less. 

The insolation of the catenary-tent-collector for 0=22. 3°N, L s =200° and albedo al=0.22 
is shown in Table 1 based on radiation data at VLI [6]. The table shows the beam, diffuse 
and the albedo insolations on sides A and B in kWhr-day and kWhr/m 2 -day. As expected, 
the beam insolation on side B is lower by 59.6% than on side A. It is interesting to note 
that the diffuse insolation comprises 46.6% of the global insolation. This characteristic is 
typical for Mars, a place where the atmosphere consists mainly of dust particles. 

Table 1. - Insolation on a catenary-tent collector, 0=22. 3°N, L s =200°, al=0.22 



BEAM 

COMPONENT 

DIFFUSE 

ALBEDO 

GLOBAL 


side A 

side B 

Total 

COMPONENT 

COMPONENT 

INSOLATION 

Q(kWhr-day) 

12.65 

5.11 

17.76 

16.27 

0.89 

34.92 

q(kWhr/m 2 -day) 

2.246 

0.907 

1.577 

. 

1.445 

0.079 

3.101 
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5. CONCLUSIONS 


The article analyses the performance of a catenary-tent-collector [1] (a flexible 
blanket that falls freely on both sides of a central support). This kind of collectors has 
characteristics (portability and simplicity) that are desirable for solar power plants on outer 
space planets. Because of its shape, there Is a self -shading effect that must be taken into 
account in the solar radiation calculation. Therefore, the shape and area of the shadow on 
the collector is calculated and used in the determination of the beam radiation. The diffuse 
and albedo radiation were also calculated to determine the global radiation on the collector. 
The numerical example is based on solar radiation data on Mars. 

6. NOMENCLATURE 


A 

Collector area, [m 2 ] 


A sh 

shaded area of a collector, [m 2 ] 


A ush 

unshaded area of a collector, [m] 


D 

length of the collector, [m] 


F A-C 

view factor of area A with respect 

to ground 

F A-S 

view factor of area A with respect 

to sky 

FdA-S 

view factor of incremental area dA 

with respect to sky 

G b 

direct beam irradiance, [W/m 2 ] 


G dh 

diffuse irradiance on a “horizontal surface, [W/m 2 ] 

G h 

global irradiance on a horizontal surface, [W/m 2 ] 

H 

height of collector, [m] 


k 

catenary constant 
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L 


collector width, [m] 


L s areocentric longitude of the sun (for Mars) 

P a/ albedo irradiance, [W] 

P* direct beam irradiance on an unshaded side of a collector, [W] 

P|* direct beam irradiance on a partially shaded side of a collector, [W] 

direct beam irradiance on an unshaded side A and B of a collector, 
respectively, [W] 

P«J,P£k direct beam irradiance on a partially shaded side A and B of a collector, 
respectively, [W] 

Px»Pj/ x and y components of a shadow length, [m] 

Q insolation, [kWhr-day] 

q insolation, [kWhr/m 2 -day] 

a sun altitude 

~t c collector azimuth 

~i s sun azimuth 

nr sr sun azimuth at sunrise 

6 solar declination angle 

e characteristic angle of the collector 

0 angle between solar ray and the normal to the collector 

4 relative shaded area 

0 local latitude 

w solar hour angle 


17 


7. REFERENCES 


[1] A.J. Colozza, "Design, Optimization, and Analysis of a Self-Deploying PV Tent 
Array”, NASA Contract Report 187119, June 1991. 

{2] F.B. Beer, E.R. Johnston, Mechanical for Engineers: Static and Dynamics 
McGraw-Hill, New YOrk, 1967, pp. 258-273. 

[3] J. Appelbaum, J. Bany, "Shadow Effect of Adjacent Solar Collectors in Large Scale 
Systems" Solar Energy, Vol. 23, pp. 497-507, 1979. 

[4] J. Bany, J. Appelbaum "The Effect of Shadowing on the Design of Field of Solar 
Collectors" Solar Cell, Vol. 20, pp. 201-228, 1987. 

[5] H.C. Hottel, A.F. Sarofim, Radiative Transfer, McGraw Hill, New York, 1967, pp. 
25-39. 

[6] J. Appelbaum D.J. Flood, "Solar Radiation on Mars", Solar Energy, Vol. 45, pp. 
353-363, 1990. 


18 


APPENDIX A 

Calculation of points E and T 


Calculation of Point E and T for case (i) 

Point E 

Point E may be calculated from the projection of the solar ray NP on the z-y plane, 
Fig. A1 (a). The projection line equation is: 

z = H(1 - y/Py) (A.l) 

The intersection of the projection line, eq. (Al) with the catenary equation f c (y) yields the 
coordinate y^, i.e., the solution of: 



The coordinate is on line OP, i.e., 

P* 

xr = rr yr 

r y 

and finally the coordinate z^ is calculated from eq. (Al). 


(A. 3) 


Point T 

Point T is determined by a sun ray that its projection on the z-y plane is tangent to 
the arc N N n , Pig- Al (b), i.e., at this point, the derivative of f c (y) is equal to the slope 
of line AB. From eq. (1) 

fc'(y) = -sinh(5_I_yj (A. 4) 
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Since AB and N P y are parallel lines, the slope of line AB is -H/P y , thus obtaining: 


y r = D - k sinh-KH/Pj,) 


(A.S) 


The coordinate xr = 0 and the coordinate z t is calculated from the f c (y) equation. 


CalculatioTi o/ yy and Xr /or case (it) 


A sun ray that touches the foot of the tent at point E passes through point U on the 
blanket, Figs. 5 and A.l.(a). It follows that 


(D - y t ,)/x c = P y /Pi 


(A. 6) 


and 


Zy _ [ciyu) _ _H 

D - yv D - y v Py 


(A. 7) 


from which y y is determined by solving: 


fc(yt') 



DH 

Py 


(A.S) 


where f c (y) is given in eq. (1). The coordinate Xe is then calculated from eq. (A. 6). 


Calculation o/ yr [or case (m ) 


Referring to Figs. 6 and A.l (b), we may write 


(y e - yi/)/L = P y /Pz 


(A. 9) 


and 


zy ~ z e _ fc(ytO ~ fcCyg) _ 
yr - Yu yr - y v Py 


(A. 10) 
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From eqs. (A. 9) and (A. 10), is determined by: 


fc(yr - l - fc(yc) = ^ 

where F c (y) is given in eq. (1). 


(A. 11) 


APPENDIX B 


Using the parabolic approximation for f c (y) greatly simplifies the mathematical 
calculation. The error is small and the results are very similar to the catenary case. The 
procedure for calculating the points for the shadows are the same as in Appendix A but now 
fc(y) is the parabolic approximation (eq. (2)), i.e.. 


f c ( y) 



. fc’(y) = ( l 



(B.l) 


Case (?) 


The coordinates of points E and T are: 



(B.2) 


(B.3) 


rs /-n 

With these results, the arcs N N n , N E y and T E y can now be determined. The length of an 
arc is given by eq. (16) and for the parabolic approximation we have: 


g(y) = j \/l + [fc’(y)l 2 dy = J yjl + ^ (D - y)j 2 


dy 


the solution of the integral is: 


6(y) = ^ sinh 1 


2H , 

^ (y - d) 


^(y - D) \[l + [j$ (y - D)] 2 (B.4) 


The arc lengths are therefore: 
iTn* = g(D) - g(0) 
N~Ei, = g(y c) - g(0) 
T Ey = g(y £ ) - g(y r ) 


(B.5) 

(B.6) 

(B.7) 


Case (») 


Solving eq. (A. 8) with eq. (2) for yv (Fig. 5) results in: 




(B.8) 


and using eq. (A. 6) we obtain the coordinates of xe- The coordinates of points E and T are: 


xr = P*(jr) 2 • yf = D , 2 e = 0 


(B.9) 


x r = 0 . y r = d(i - , z r = h(^) 2 


(B. 10) 


The expression for T N„ is given by; 


T N„ = B(yr) " 6(D) 


(B.ll) 
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where g(y) is given in eq. (4). 


Case (m) 

Solving eq. (A. 11) with eq. (2) for (Fig. 6) results in: 

^ = D ( 1_ ft) + L fe (B12) 

The coordinates of points E and T are: 
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Table 1. - Insolation on a catenary-tent collector, 0 = 22. 3°N, L s =200°, al=0.22 



BEAM 

COMPONENT 

DIFFUSE 

ALBEDO 

GLOBAL 


side A 

side B 

Total 

COMPONENT 

COMPONENT 

INSOLATION 

Q(kWhr-day) 

12.65 

5.11 

17.76 

16.27 

0.89 

34.92 

q(kWhr/m 2 -day) 

2.246 

0.907 

1.577 

1.445 

0.079 

3.101 



Figure 1 . — Artist's view of a catenary-tent-collector. 
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Figure 8. — Catenary collector - solar radiation calculation on an unshaded side. 
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